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Abstract 

Motivated by recent transport and neutron scattering experiments suggesting an orientational symmetry 
breaking in underdoped cuprates we present a theoretical study of Raman scattering near a cf-wave Pomer- 
anchuk instability (PI). The J-wave component of Raman scattering from electrons and phonons allows to 
study directly order parameter fluctuations associated with the PI. Approaching the PI from the normal state 
by lowering the temperature a central peak emerges both in electronic and, as an additional low-frequency 
feature, in phononic scattering. Approaching the PI in the superconducting state at low temperature by 
decreasing the doping concentration the central peak is replaced by a soft mode with strongly decreasing 
width and energy and increasing spectral weight. These predicted low-energy features in Raman scattering 
could confirm in a rather direct way the presence of a PI in high-temperature cuprate superconductors and 
in Sr3Ru207. 

PACS numbers: 74.25.nd, 74.25.Kc, 71.18.+y, 74.72.-h 
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I. INTRODUCTION 



In condensed matter electrons move through a crystal lattice whose symmetry is characterized 
by a point group. The electronic band structure usually has the same symmetry as the lattice 
and so does the Fermi surface. However, it was shown that the symmetry of the Fermi surface 
can be broken spontaneously by electron-electron correlations in the two-dimensional t-J-- and 
Hubbard 4 5 models leading to spontaneous Fermi surface deformations characterized by a J-wave 
symmetry [rf-wave Fermi surface deformations (JFSD)]. This instability is frequently referred 
to as a rf-wave Pomeranchuk instability, which is characterized by the violation of the stability 
criteria for isotropic Fermi liquids derived by Pomeranchuk. 6 However, it should be noted that 
the JFSD state can also be realized not only in strongly correlated electron systems such as those 
described by the t-J model 1 3 but also without a breaking of Pomeranchuk's stability criterion 
in systems where the transition can become of first order at low temperatures. 7 8 The JFSD state 
breaks only the orientational symmetry, that is, its instability is driven by zero momentum charge- 
density fluctuations with internal <i-wave symmetry and leads to an electronic nematic state. As 
originally introduced in Ref. |^, an electronic nematic state can also be realized by invoking charge 
stripes. In a first step both orientational and translational symmetry are broken by condensing the 
electrons into a charge stripe state characterized by a set of large wave vectors which break the 
orientational symmetry. In a second step the stripes melt restoring the translational but not the 
orientational symmetry. In the following we restrict ourselves to the case where the JFSD leads 
directly to an electronic nematic state without showing first an instability towards stripes. 

The double-layer strontium ruthenate Sr 3 Ru 2 7 (Sr327) has attracted much attention as a 
compound likely exhibiting a JFSD instability. 1 ^— Compelling, but indirect, evidence for this 
comes from the observation of a strong xy anisotropy of the resistivity which is present only 
in the ordered phase. 11 Angle-resolved photoemission spectroscopy (ARPES) 13 and de Hass-van 
Alphen 1416 measurements could detect Fermi surface deformations directly, but convincing ex- 
perimental evidence for their existence has not been obtained yet. Theoretically many properties 
have been successfully interpreted in terms of a JFSD instability, for instance, the metamagnetic 
transition, 17 the enhancement of the residual resistivity, 18 the phase diagram and various thermody- 
namic quantities, 19 universal numbers, 20 the bilayer effect, 21,22 suppression of a critical temperature 
due to impurities, 23 the spin-orbit effect, 24 and orbital degree of freedom. 25 Theoretical predictions 
based on the JFSD instability were also made for the pattern of Friedel oscillation around an 
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impurity,— the attenuation of ultrasound waves, 27 and the singular behavior of the uniform mag- 
netic susceptibility at the JFSD instability.— 

In the case of the high-temperature superconductors YBa 2 Cu 3 O v (YBCO v ), the dynamical in- 
plane magnetic susceptibility is strongly anisotropic, both for slightly underdoped (YBCC^.e)^ - 
and optimally doped (YBCOg.ss)^ compounds. The anisotropy increases with decreasing doping 
and is most pronounced around the onset temperature of superconductivity or of the pseudogap, 
whereas it is suppressed in the superconducting state. It was shown theoretically that these fea- 
tures can be well understood in terms of the competition of the singlet pairing formation and JFSD 
correlations.— In the strongly underdoped region (YBCO6.45) neutron scattering experiments re- 
vealed a qualitatively different feature of the anisotropy. 32 The in-plane anisotropy of the magnetic 
excitation spectrum increases monotonically below 150 K, saturates below 50 K, but is not sup- 
pressed below T c = 35 K. Moreover, the low-energy spectral weight does not decrease below 
T c but is rather enhanced. These peculiar phenomena can be interpreted as i) a quantum phase 
transition to the JFSD state deep inside the superconducting state^i or ii) a substantial suppres- 
sion of singlet pairing due to the competition with increasing JFSD correlations in the strongly 
underdoped region.— 

Quite recently the measurement of the Nernst coefficient in the doping region from 1 1 - 18 % 
in YBCO^ 6 . showed a strong xy anisotropy. It sets in near the temperature where the pseudogap 
appears so that the pseudogap region is interpreted as the region with a finite JFSD in agreement 
with a theoretical study. 37 However, one should note that the regions where the in-plane anisotropy 
has been observed by neutron scattering 29130 ' 32 and by transport 36 differ from each other so that it 
is difficult at present to reach clear-cut conclusions. The experimental evidence for nematic order 
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in cuprates has recently been critically reviewed in Ref. 

Usually an emergent instability can be studied by measuring the enhancement of the 
corresponding susceptibility. The susceptibility describing the JFSD is the J-wave charge 
compressibility,*^ which can be measured directly by Raman scattering. Hence Raman scat- 
tering can provide decisive evidence for a JFSD instability and its correlations in actual systems. 
However, despite various experimental studies in Sr327 and YBCO, Raman scattering experiments 
have not been reported to confirm a JFSD in those materials. 

In this paper we provide theoretical predictions of the Raman scattering intensity from electrons 
and phonons near the JFSD instability in both normal and superconducting states by employing 
parameters appropriate to cuprate superconductors. In the superconducting state the Raman scat- 
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tering intensity can be computed in terms of the non-interacting electron propagator, i.e., without 
considering the damping of electrons. In the normal state, however, it is crucially important to 
include the electronic self-energy. We therefore include the Fock diagram for the self-energy, 
express it in terms of the bosonic spectral function a 2 F(co) and fit the latter to the self-energy 
measured in ARPES.^i° 

The paper is structured as follows. In Sec. II we present formulas for electronic and phononic 
Raman scattering near a JFSD instability. Since the order parameter for JFSD fluctuations has 
Big symmetry for a square lattice only the Bi g component of the Raman tensor and Bi g phonons 
will be directly affected by order parameter fluctuations. In Sec. Ill we study Raman scattering for 
two different ways to approach the JFSD instability. In the first case the system is always in the 
normal state and the temperature is lowered for a fixed doping in the underdoped region. In the 
second case we assume that at low temperature the JFSD instability lies in the superconducting 
state and is reached by decreasing the doping. Results for both cases are given in this section. 
Sec. IV contains a detailed discussion of these results and our conclusions. 



II. FORMALISM 



In the following we will consider fermions on a square lattice which has the tetragonal point 
group symmetry D 4h . Since the order parameter of the JFSD is the charge density with internal 
J-wave symmetry and zero total momentum, JFSD fluctuations will be most easily detected in 
the B lg component of Raman scattering and for a zone-center phonon with B lg symmetry. We 
therefore will focus on these two quantities in the following. Throughout the paper we will also 
use the lattice constant of the square lattice as the length unit. 



A. Electronic Raman scattering 

The electronic contribution to the B lg Raman vertex is given in the effective mass 
approximation^ 

7k 2\dk 2 x d**)' K) 
where e k is the electronic dispersion, 

6 k = -2?(cos k x + cos k y ) - At' cos k x cos k y - 2?"(cos 2k x + cos 2k y ) , (2) 
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with t, t' , and t" being the nearest, second-nearest, and third-nearest neighbor hopping integrals, 
respectively. Inserting Eq. © into Eq. (Q~|) yields 

y B>g = ?(cos k x - cos k y ) [l + St" (cos k x + cos k y )/t\ . (3) 

The Raman scattering intensity S (oS) is given by 

S(fo) = --[1 + b(co)]Im X B Hu) , (4) 
n 

where b(oS) is the Bose function given by (e^ - IT 1 and = T is the temperature. The quantity 
;f Blg (a>) is the retarded Green's function with two Raman vertices as end points and is given by 

^(co) = ~ dtc^ +mt ([p d (t),p d (0)]), (5) 
Jo 

where /V is the total number of lattice sites, + an infinitesimally small quantity; <• • •) denotes the 
equilibrium expectation value, [•, •] is the commutator, and p d (t) is the Heisenberg representation 
of the J-wave charge density operator 

Pd = J] rl >gc L c ^ (6) 

k,o- 

with c^ cr (cko-) being the creation (annihilation) operator of electrons with spin cr and momentum 
k. Within the RPA x Bls (w) is given by the bubble diagrams shown in Figs. CUa) and (b). The 
double line represents the electronic Green's function. In the normal state, which we discuss first, 
self-energy corrections in the electronic Green's functions must be taken into account, as shown 
diagrammatically in Fig.QJc). Otherwise each bubble would become zero in the zero momentum 
limit at every finite frequency. This means that mainly the incoherent part of the Green's function 
contributes to Raman scattering in the normal state. 

In order to get a finite self-energy we consider the coupling of electrons to some bosonic fluc- 
tuations, described by the Fock diagram shown in Fig. [Del). Analytically one obtains, adopting 
the usual approximations in evaluating the Eliashberg equations, 42 

ImSM = -7T I dv a 2 F(v) [2b(v) + f(v - to) + f(v + io)] , (7) 
Jo 

where f(v) = (e 6 " + l)" 1 is the Fermi function, and a 2 F(v) specifies the bosonic spectral func- 
tion. We have neglected the momentum dependence of a 2 F for simplicity. Note that although 
the notation of a 2 F(v) is often used in the context of a phonon spectrum, the bosonic modes in 
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(a) X Bl8 =<^> + «^^^^^ 




(c) = = + ~ @H= 

(d) © = £% 

FIG. 1: (a) Graphical representation of x Big . The vertex with a circle (square) indicates the form factor 
j\ lg (dk)- (b) Effective electron-electron interaction driving the <iFSD instability, (c) Full electronic Green's 
function. The single solid line denotes the free electron propagator with the dispersion e^. (d) Electronic 
self-energy originating from the coupling to some bosonic fluctuations represented by the sawlike line. 



Fig. Old) are arbitrary in our model. We model the function a 2 F(v) with three parameters, a , vo, 
and v c , as illustrated in Fig. |2ta). Comparing with ARPES measurements in the normal state in 
cuprates 39 40 we choose v = v c /4, v c = 2t/3, and a = 1/4 with t « 150 meV. In Fig. |2tb) we 
show lmL(oS) for several choices of T. The magnitude of lmL(to) becomes larger with increasing 
T, indicating reductions of the lifetime of quasi-particles at higher T. As a function of energy, 
on the other hand, the longest lifetime of quasi-particle is realized on the Fermi surface, namely 
at to = 0; the magnitude of Im2(o;) is enhanced with increasing to and saturates to the value, 
-7i dv a 2 F(v)[2b(v) + 1] for to — > oo. The real part of is computed numerically from the 
Kramers-Kronig relation, Re2(o;) = jP-V. dv I ^ v) , where the integral is defined as the prin- 
ciple value denoted by "P.V.". The obtained ReE(a») is shown in Fig. 0c). The real part of 
vanishes at to = 0. Its magnitude forms a peak around to « 0.5t and is suppressed at high to with 
a tail characterized by to' 1 . Our self-energy reproduces well the data 3 - 9 -^ extracted from ARPES 
measurements in cuprate superconductors. The spectral function of the full Green's function 
[Fig. djc)] is given by 

A(k,to) = — ImG(k, to) (8) 

_ 1 (9) 

n[to- (£ k - n) - ReE(a))] 2 + [ImZ(w)] 2 ' 

Here fi is the chemical potential which is approximately determined by the relation 6=1 — 

ji Zk f( e k ~ H) f° r a given doping concentration 5 and T. Figure [3£a) shows A(k, to) as a function 

of to for several momenta along the (0, 0) - (n, n) direction. A relatively sharp peak is seen only 
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FIG. 2: (a) Model of a 2 F(v). (b) Imaginary part and (c) real part of the electronic self-energy for several 
choices of T in the normal state for vq - v c /4, v c = 2/3, and ao = 1 /4. The energy unit is taken as t. 




FIG. 3: (a) A(k, oS) as a function of oj at the momentum k - = 2nk r {\, 1) with k r ranging from zero 
to -0.10 with an interval of 0.01; kp is the Fermi momentum along the (0,0) - (n, n) direction, (b) The 
renormalized electronic dispersion. 



close to the Fermi energy (u> = 0), away from the Fermi surface it is substantially broadened 
because of the presence of the sizable ImL(co). The peak position of A(k, oS) is plotted in Fig. Ob). 
The renormalized electronic band dispersion displays a kink, as indicated by an arrow, due to the 
coupling to the bosonic fluctuations [Fig.[TJd)]. These features in Fig.|3]are qualitatively consistent 
with ARPES data.^i2 
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A single bubble diagram in Fig. [2 a) corresponds to the analytical expression, 

ITV) = Tr Z f deide 2 A(k, ei )A(k, e 2 ) } ~ /fe) n+ , (10) 
N ^ J-oo 61 + oj - 6 2 + i0 + 

where the form factors of the vertices are denoted by a k and /? k . For the imaginary part lmll al3 {a)) 
we obtain 

ImIP"(w) = ^ ttw 8 k I deA(k, e)A(k, 6 + w) [/(e + <u) - /(e)] . (11) 

^ %J — 00 

The real part Ren Q ^(o>) is determined from the Kramers-Kronig relation 

Ren^) = I P .v. fdv555!!W. (12) 

^ J-co v - a» 
Finally the Raman response function ^f Blg (o;), described in Fig. HJa), is given by 

x b Hoj) = u^(oj) + IFW — ^--n^M . (13) 

1 - gU da {a)) 

The superscripts of II(a») "7" and "d" indicate the form factors of the vertices of the bubble dia- 
gram, which are taken as y^_ lg [Eq. ©] and d k = cos k x - cos k y , respectively. The rf-wave form 
factor comes from the electron-electron interaction which we write as 

^ E 8 d k d * c L- c *°- c l'o-' c v°-' ' ( 14 ) 

kk' era' 

where g(< 0) is the coupling strength. This interaction generates the effective interaction shown 
in Fig. \V[b) and drives the JFSD instability, as was extensively studied theoretically . 1 i 4 i 7 8 i 43 The 
condition for the instability is given by 

l-gIT M (0) = 0. (15) 

From Eqs. © and (fTTI)- (fT3l we computed the Raman scattering intensity numerically employing 
the self-energy shown in Fig. |2] 

The selection of diagrams in Fig. \T\ corresponds to the lowest-order conserving approximation 
in the sense of Baym and Kadanoff. 44 The diagrams shown in Fig. 3 in this reference also apply 
in our case if we consider the dashed line as a sum of the interaction of our Eq. (TT4)) and our 
boson-mediated, retarded interaction which we have assumed to be independent of momentum. 
The Hartree terms to the self-energy can be omitted because they are either zero or represent 
just a renormalization of the chemical potential. The interaction Eq. ([141) does not contribute in 
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the thermodynamic limit to the Fock term of the self-energy in contrast to the boson-mediated 
interaction which yielded the contribution given in Eq. ©. The vertex is given as the functional 
derivative of the self-energy with respect to the Green's function. Limiting ourselves to the J- wave 
vertex we see that it is only the functional derivative of the Hartree term of the interaction Eq. (fT4l) 
which contributes to the vertex and produces the chain of bubbles in the rf-wave susceptibility. 
This means that our approximation scheme respects all conservation laws and should be free of 
artifacts due to an inconsistent approximation. 

In the superconducting state the quasi-particle contribution to the Raman scattering intensity is 
finite at finite frequencies. Since the self-energies are also much smaller in the superconducting 
state compared to those in the normal state it seems to be reasonable to neglect self-energy effects 
in this case. 41 Assuming the following form for the J- wave superconducting gap 

A k = ^A (cos k x - cos k y ) (16) 

and the band dispersion Eq. ©, we obtain for the single bubble diagram [Fig.QJa) with S = 0], 

irV) = i X ^4 tanh ( * T - * +T ) . (iv) 

N ^ 2 \(0- 2E k + iY co + 2E k + iT/ 

where E k = ^(e k - jj) 2 + A^. Y is an infinitesimally small positive quantity which we approximate 
in our numerical calculations by Y = 0.001?. The chemical potential is approximately determined 
from the relation 6 = jj 2k tanh The Raman scattering intensity and response function 
are given again by the formulas Eqs. © and (fT3l) . respectively. 



B. Raman scattering from phonons 

Raman scattering can also determine the spectral function of phonons. Since the interaction 
driving the JFSD instability couples to phonons with B ig symmetry, we focus on phonons with 
this symmetry. The corresponding electron-phonon coupling contains the J-wave form factor d k 
and is given by 

gkk = g P hd k , (18) 

where g kk is the electron-phonon matrix element for the electronic momentum k and vanishing 
momentum for the phonon; g ph is the coupling constant. The non-interacting retarded phonon 
propagator is given by 

1 1 




FIG. 4: Graphical representation of the phonon propagator (double dashed line); the single dashed line 
denotes the non-interacting phonon propagator, the solid square the form factor g p hdk- The rest of the 
notation is the same as Fig. [TJ 



where co is the energy of the zero-momentum B lg phonon, which corresponds to co = 4f/15(~ 40 
meV for t » 150 meV) in YBCO.— The full phonon propagator is given graphically in Fig.BJa), 
namely, 

D- 1 M = D- 1 M-£ ph M. (20) 

The free phonon propagator is renormalized by the electron-phonon interaction, which picks up 
the correlation function of the JFSD instability as shown in Fig. @Jb). The phonon self-energy 
2ph(<^) [Fig.BJb)] has exactly the same structure as Fig. [2 a) except for the difference of vertices. 
The computation of S ph (a») is straightforward, yielding 

W =s, H- g u^) (21) 

= gl h U dd (co), (22) 



where we have introduced 



(23) 



The quantity Ti dd (a>) would become identical to x B[g (of) if "y" were replaced by "J" in Eq. (TOT) 
The Raman intensity 5 P h for phonon scattering becomes 

1 



Sph(w) = — [1 + b(io)]lmD(to) , (24) 
n 



where from Eqs. (|20l ) and (1221) . 



! t nr , 40J ° 8 ^ lmU dd (oj) 

-lmD(to) = — — . (25) 

* * [co 2 - co 2 - 2co gl h RcU dd (co)] 2 + [2co () gl h lmU dd (co)] 2 

Since 11^(0*) has already been computed both in the normal and superconducting state in Eqs. (fTTI) . 
(fT2l) . and (fTTT) . the Raman intensity S ph (co) is easily obtained from Eqs. (f23T)-(f25l). 
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1. Renormalization of the dFSD by the electron-phonon coupling 



It is instructive to provide an expression of the static rf-wave charge compressibility M^- the 
susceptibility associated with the dFSD instability. This quantity is obtained by summing up 
the bubble diagrams connected by electron-electron and electron-phonon interactions as shown 
graphically in Fig. [51 that is, 

n**(0) n ~ 

K d = i t, — (26) 

l-[g + g2 h D (0)]n^(0) 

n^(0) 

(27) 



i -gu dd (oy 

where 

2^ 2 



0)Q 

,2 



(28) 



is a renormalized coupling constant. Since both g^ h and a>o are positive the original interaction 
g(< 0) is enhanced to become \g\ > \g\. The coupling to the B ig phonon mode therefore increases 
the attractive interaction causing the dFSD instability by the amount — . Therefore the dFSD 



instability can occur more easily and the instability condition Eq. (1151) is replaced by 

1 - gU dd (0) = . (29) 




-^jyv = -wv\, + 

FIG. 5: Graphical representation of the <i-wave charge compressibility. A spring denotes two interactions, 
the electron-electron interaction (wavy line) and the electron-phonon interaction (dashed line), and the 
corresponding form factors at the vertex (open circle) are and g p h<4, respectively. The shaded vertex is 
defined by the equality of the first and second line, the rest of the notation is the same as in Fig. [TJ 



11 



2. Strength of the electron-phonon coupling 



An estimate for the coupling constant g ph is obtained both from first-principle calculations 
based on the local density approximation (LDA) and from experiment. In general, the dimension- 
less coupling constant A for a phonon with energy co and zero momentum is defined by 

A = 2jV(0)<|£ kk | 2 > FS M, (30) 

where N(0) is the density of states at the Fermi energy for one spin direction, (• ■ -) FS denotes an 
average over the Fermi surface, and g kk is defined by Eq. (TTST ). First-principle LDA calculations 
yielded for the Bi g phonon mode with 40 meV the values A = 0.02 46 and 0.06. 47 On the other hand 
A is related to superconductivity-induced self-energy effects of the phonon. It has been argued that 
the observed self-energy effects are compatible with these values for A, especially, with the first 
value. 46 A simple connection between theory and experiment can be also obtained by noting that a 
phonon well below the superconducting gap at T = should show a relative frequency softening of 
about Stoo/coo = A. 46 The above 40 meV phonon softens by ~ 1 meV due to superconductivity 4 -^ 
yielding A « 0.02 - 0.03 in rough agreement with the theoretical prediction. 

For a phonon one has gkk = g P h<4 and (d 2 )^ = Nd(0)/N(0) where Na(0) is the d- wave 
projected density at the Fermi energy, namely A^(0) = Jd 2 kc>(ek - pi)d^/(2n) 2 . We thus obtain 

Acoq 

gph 2N d (0) K } 

Since 2A^(0) is equal to the low-energy limit of a single bubble at T = 0, namely -IT^O), we find 
2N d (0) ~ l/\g\, so that g 2 h ~ OAAt 2 for co = At/ 15 and g = -1.5?, yielding values between 0.008 
and 0.024 for g 2 h . In our numerical calculations we use the representative value 0.02. 



HI. RESULTS 

Guided by experiments in YBCO^ 6 - we would like to choose one parameter set in our simple 
model such that the JFSD instability is reached (a) with decreasing temperature at around T ~ t/10 
in the normal state with a doping concentration 6 = 0.10 and (b) with decreasing doping at around 
6 ~ 0.20 for r«0. These conditions are approximately fulfilled for f = t" = and g = -1.5? in 
our model. 50 For convenience, we use t as the energy unit in presenting our results. Experimentally, 
the effective t has a value of about 150 meV. 
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A. Electronic Raman scattering 



In the normal state we fix the doping to 6 = 0.10 and consider the temperature as a tuning 
parameter to approach the JFSD instability from high temperatures. For our parameters the JFSD 
instability occurs at T = 0.098. In Fig. [(J a) we show the co dependence of lmx Blg (co) for a sequence 
of temperatures T ranging from 0.10 to 0.20. At high T the weight of Im^ Blg (o>) extends very 
broadly over the whole energy region shown in Fig. [(J a). With decreasing T the low energy weight 
(co < 0.2) gradually increases and sharpens up to form a very steep peak near zero frequency. In 
Fig. (6^b) we plot the function S(ca), defined in Eq. ©, which is measured in a Raman scattering 
experiment. Although the peak position is not exactly at co = 0, S (co) displays essentially a central 
peak already well away from the critical temperature. Its spectral weight increases strongly when 
the critical temperature is approached from above. The energy dependence of lmn dd (co) is shown 
in Fig. (6fc) for several values of T. While lmn dd (co) exhibits also a pronounced peak its energy 
is much larger than that of Im^ Bl «(a»). Moreover, the effect of temperature is much weaker in 
lmn dd (co) than in Im^ Blg (&>). The real part of Tl dd (co) is shown in Fig.[6td). Its magnitude forms a 
broad peak at co = at high temperatures which sharpens up with decreasing temperature. Since 
the JFSD instability occurs when Eq. (fT5l) is fulfilled, collective fluctuations of the JFSD develop 
when the magnitude of ReU dd (0) approaches l/|g| = 2/3 with decreasing T. Hence the very 
pronounced peak of lmx B[g (co) at low energy, seen in Fig. [(J a X is a direct consequence of the 
development of JFSD correlations. 

Next we investigate the evolution of JFSD fluctuations in the superconducting state at T = 
0.001. Here we take the doping concentration 6 as a tuning parameter and approach the JFSD 
instability by decreasing 6. We choose the superconducting gap amplitude to be A = t/5 which 
seems to be reasonable for cuprate superconductors; the doping dependence of A is not important 
for our conclusions and will not be considered. We obtain 6 C = 0.207 for the critical doping rate 
where the JFSD instability occurs and consider the region 6 > 6 C . Since there is no qualitative 
difference between S(co) and Im^ Blg (o>) we present results only for S(co). In Fig. Ha) we show 
S(co) as a function of co for several choices of doping; similarly Imn^(a>) and ReU dd (co) are 
presented in Fig.Hb) at 6 = 0.25 and 0.50. At 6 = 0.50 the peak position of S(co) is nearly the 
same as that of ImU dd (co) because the peak originates from individual excitations. Its position 
is determined approximately by co = 2|A k | at k = k F = (kp, 0) or (0, kp), where k F is the Fermi 
momentum along the k x or k y direction. With decreasing 6 the peak of S (co) shifts to lower energies, 
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FIG. 6: (Color online) (a) to dependence of lmx Bl °(a>) for a sequence of temperatures T close to the cfFSD 
instability at T c = 0.098 in the normal state at 8 = 0.10. (b) Low energy region of S(to) near the <fFSD 
instability, (c) lmH dd (to) and (d) ReIl dd (to) as a function of to for several values of T. 



its half-width decreases and its height increases strongly. The peak position substantially deviates 
from that of ImIT w (a>) [see the results at 6 = 0.25 in Figs[7|a) and (b)], indicating the development 
of collective fluctuations of the JFSD. In fact, the peak position of S(a>) is near the instability 
determined by the resonance condition, 

1 - sRelT'Vres) = . (32) 

The resonance energy a> Tes is plotted in Fig.[7£c) together with the peak energy of S (a») and its peak 
height. At high 6, Eq. ([32b does not have a solution and the peak of S (of) must be attributed mainly 
to individual excitations. For 6 < 0.40 Eq. (1321) has a solution. It is seen that upon approaching 
S c , k> P eak becomes almost identical with o» res . Since lmn dd {oS) « at a> « o> res the evolution of 
S((o) in Fig.|7]^a) indicates the development of a well-defined collective mode associated with the 
JFSD. Because of the collective fluctuations the peak intensity of S (to) is strongly enhanced upon 
approaching 6 C and diverges at 6 = 8 C . The peak energy vanishes as a) peak ~ (6 - 8 C ) 112 , which can 
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FIG. 7: (Color online) (a) co dependence of S(co) for a sequence of doping concentrations in the supercon- 
ducting state at T = 0.001; the actual value of S(a>) is obtained by multiplication with the factor indicated 
near each peak except for 8 = 0.50. (b) co dependence of ReYl dd (to) (solid line) and lmH dd (co) (dashed 
line) at 6 = 0.25 and 0.50. (c) The peak position of S(co) (solid circles) and its peak height (solid line) 
as a function of doping; also shown are the energies a> res (open circles); the <iFSD instability occurs at 
5 C = 0.207. 



be read off from Fig. |7£c). 

It is interesting to note the different evolution of S (co) in the normal and the superconducting 
state. In the normal state the magnitude of ReH dd (a>) has a maximum at co = and decreases 
with co [Fig. [(Jd)], whereas in the superconducting state the magnitude of RelT^O) corresponds 
to a local minimum and increases with co [Fig. |7£b)]. In contrast to the superconducting case the 
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resonance condition Eq. (1321) is not fulfilled in the normal state except at T = T c and to = 0. 
This explains why S (to) develops a central peak in the normal [Fig. Oa)] and a soft mode in the 
superconducting [Fig. Ha)] state and why the width of the peaks is much larger in the normal than 
in the superconducting state. 

B. Raman scattering from phonons 

Raman scattering from B lg phonons exhibits characteristic features near the JFSD instability. 
As a prominent example we consider the 40 meV phonon in YBCO 45 which has in an approximate 
tetragonal classification, where the chains are neglected, B ig symmetry. Our parameter values 
become to = 4/15 and g 2 h = 0.02, as discussed in Sec. II.B.2. 

Figure Oa) shows S P h(w) in the normal state for 6 = 0.10 and several values for T. In the 
presence of the electron-phonon interaction the critical temperature occurs at T c = 0.126 which 
is higher than in the case without electron-phonon interaction. Well above this temperature, for 
instance, at T = 0.20 S ph (to) consists of one single peak at to ~ cjq representing a usual quasi- 
harmonic phonon. The peak position moves only slightly upwards, its spectral weight decreases 
somewhat with decreasing temperature. However, at low frequencies dramatic changes occur: 
approaching T c from high temperatures a central peak develops. It extends over a rather broad 
energy region co < 0.1, but nearer to the instability its half-width decreases and its spectral weight 
increases strongly. It is caused by the coupling of the phonon to JFSD fluctuations. The occurrence 
of a double-peak in the phonon spectral function can be understood by studying the denominator 
of the phonon spectral function Eq. (1231) . omitting lmIl dd (co), 

Res(w) = co 2 - co 2 ~ 2co g 2 ph Rdn dd (co) . (33) 

As shown in Fig. E£b), Re.H dd (to) decreases monotonically with frequency and becomes small 
around to « ojq. Because g 2 h = 0.02 is also small Res(&>) becomes zero at co ~ co giving rise 
to the quasi-harmonic phonon mode. Since the magnitude of ReIl dd (co) assumes its maximum at 
co = and increases there with decreasing T, it eventually reaches the value co /(2g 2 h ) so that 
Res(a») becomes zero also at to = 0. This situation occurs just at T c , because the expression 
ReH dd (0) = -too/(2g 2 h ) reduces to RelT^O) = 1/g via Eq. (1231) which corresponds to the onset 
of the JFSD instability [Eq. (|29l)l. Hence both Res(a») and ImT[ dd (to) for co ~ are very small 
near T c which causes a central peak close to the JFSD instability. The long tail in frequency of the 
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FIG. 8: (Color online) co dependence of S p h (a) and Retl dd (oj) (b) for several choices of temperatures in the 
normal state at 8 = 0.10; the renormalized critical temperature of the JFSD instability is T c = 0.126. (c) 
Peak positions co ph and ^ h (< a> ph ) of S p h(co) and their peak heights as a function of T. 



central peak reflects the fact that the magnitudes of Res(a/) and Imll (oj) only slowly increase with 
increasing to. In Fig. [S£c) we plotted the peak positions and peak heights of the phonon spectral 
function as a function of T. The lower peak position is denoted by tof 1 . We see that the central 
peak emerges well above T c = 0.126 and acquires rapidly a large spectral weight with decreasing 
T which diverges at T = T c . The high-frequency part of the phonon spectral function does not 
show a pronounced temperature dependence despite the proximity to the JFSD instability. The 
peak intensity near co = ojq is suppressed at lower T because of the increase of the magnitude of 
ImIT irf (a>) around a> « a> upon approaching T c . 
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FIG. 9: (Color online) (a) to dependence of S p h(w) for several choices of 8 in the superconducting state at 
T = 0.001; the actual value at 6 = 0.50 is obtained by multiplying with a factor of two. (b) co dependence 
of Res(a») for different 6. (c) Peak positions a» ph and a> ph (< co ph ) of S p h{co) as a function of 6 together with 
their peak heights. The <iFSD instability occurs at 8 C = 0.233. 



The phonon spectral function in the superconducting state is shown in Fig. H^a) for several 
doping concentrations. At 5 = 0.50 the spectral function shows a quasi-harmonic phonon with 
one sharp peak at co '=. co\ « oj , where Res(a»i) = [Eq. (1331) and Fig. |9fb)] . With decreasing 
8 the position of this peak shifts only slightly to higher energies but is essentially unchanged. At 
5 = 0.30 an additional broad peak emerges at the low energy co « 0.14 in Fig. [9][a). Its position 
is approximately given by the energy at which Res(a>) forms a local maximum, see Fig. |9tb). 
Decreasing 6 further to 0.25 the lower peak becomes sharper, moves to lower energies and its 
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spectral weight increases. The equation Res(a>) = has now three solutions, co\, u>2, and toj, see 
Fig-Hb), with a>\ > co 2 > ^3 > 0. The solution coy « coq yields the sharp high-frequency peak in 
Sph(a)), while the solution co 3 is responsible for the low-frequency peak. The solution of co 2 cannot 
produce a peak in S p h(w) because lmX\ dd (oS) has a peak near the energy to 2 and thus broadens out 
any structure in this frequency region. These three solutions exist until o> 3 becomes zero. At this 
point the JFSD instability occurs which follows from a similar argument as given below Eq. (1331) . 
In Fig. |9fc) we present the peak positions of 5 P h(o>) and their heights as a function of 6. The 
high-frequency peak in the phonon spectral function, appearing around co « co = 4/15, displays 
only a very weak doping dependence in spite of the proximity to the JFSD instability. Its height 
also depends only weakly on 6 on a logarithmic scale. The coupling of the phonon to collective 
fluctuations of the JFSD leads to the appearance of a second peak at low energy which softens 
in frequency and increases in intensity upon approaching the JFSD instability. The lower peak 
energy vanishes as ~ (6 - 6 c ) l/2 and its intensity diverges when approaching the critical doping 
8 C = 0.233. The emergent low-energy peak is a well-defined collective mode driven by fluctuations 
of the JFSD in the sense that the resonance condition Res(a») « as well as lmYl dd (oj) « is 
fulfilled at the peak energy. 

It is intriguing to realize that the original quasi-harmonic B lg phonon mode does not behave like 
a soft-phonon when the JFSD instability is approached. Instead the phonon spectral function splits 
into a high-frequency part which is practically unaffected by the instability and an emergent low- 
frequency part which behaves like a soft and a central mode in the superconducting and normal 
state, respectively. This double-peak structure in the spectral function is robust if the original 
phonon energy cjq is sufficiently large. Otherwise, the phonon peak in the normal state may overlap 
with the emergent low-energy structure broadening the double-peak structure into a seemingly 
single peak. In the superconducting state a»o should be chosen to be larger than the peak energy 
in Imn <w (w), which is approximately given by the peak position of S(a>) shown in Fig. [7ja). 
Otherwise, the original phonon mode softens down to zero energy upon approaching the JFSD 
instability and no additional low-energy peak emerges, in contrast to Fig.[9£a). 

IV. DISCUSSION AND CONCLUSION 

We have studied Raman scattering in a system where the interaction between electrons drives 
the system towards a JFSD instability, both in the normal and superconducting state. The electrons 
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are assumed to live on a square lattice with hopping amplitudes t, f and t" between first, second 
and third nearest neighbors. The interaction is a charge-density interaction with internal J-wave 
symmetry and interaction strength g. The parameters t' , t" , and g are set up to mimic the strong 
tendency towards the JFSD instability in YBCO. One could wonder whether our choice of f = 
t" = and g = —l.5t is unrealistic because the presence of substantial second- and third-nearest 
neighbor hoppings is well known in cuprates 51 and the value of g simply seems too big. We would 
like to stress that the above parameter values should be interpreted as effective parameters within a 
phenomenological approach. 50 It is worth mentioning that a mean-field JFSD instability occurs in 
the t-J model with realistic parameters t', t", and g(= -37/8) for cuprates at lower carrier densities 
at temperatures as high as 0.27.J However, it is not easy to perform the above calculations directly 
for the t-J model. Nevertheless, we believe that the essential features of such a more microscopic 
approach are retained at least qualitatively in our simple phenomenological treatment. 

One result of our calculation is that the spectral function of a B lg phonon exhibits a double- 
peak structure when the JFSD instability is approached as shown in Figs. [8] and [9] The double-peak 
structure might seem similar to the emergence of a central peak near structural phase transitions 
for several perovskites such as SrTi03,— & LaA103,- 4 and KMnF 3 . 55 However the quasi-harmonic 
phonon also exhibits softening for these materials, in contrast to our results. Moreover the central 
mode in the experiments has been explained in terms of impurity scattering, 56 a different mecha- 
nism from ours. The double-peak structure we have obtained in Figs.[8]and|9]can be interpreted as 
a general aspect in a coupled system of phonons and order parameter fluctuations. In fact, similar 
results to ours were obtained in a different context, for example, in pseudospin-phonon systems 57 
and in superconductors with a strong electron-phonon coupling 58 59 explaining the double peaks 
of Raman spectrum with E 2g symmetry observed in MgB 2 .— 

For YBCO a strong tendency toward xy symmetry breaking was observed. 29 ' 30 ' 32 ' 36 Its order 
parameter may be defined by 

<P = \ J]< c Lr c * + clc i+xtT ) - (c] +y(T c i(T + clc i+y(T ) (34) 

= ^ d k {clc k(r ) , (35) 

k<r 

where i denotes the site on a square lattice and we have assumed that (f> is constant. Equation (1341) 
or (|35l) is nothing but the order parameter of a JFSD instability.-i^ 8 - It is characterized by Ising 
symmetry and thus two solutions, (p = <Po an d ~4>o, are degenerate. In order to favor either solution, 
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it may be natural to apply a small external perturbation which breaks xy symmetry in the Cu02 
plane. In fact the compound YBCO contains the CuO chains, which serves as a uniaxial strain. In 
this case, the JFSD instability becomes a crossover phenomenon, but the crossover is still sharp 
as far as the external anisotropy is weak which seems to hold in YBCO. 

From the very strong anisotropy of the magnetic excitation spectrum in YBCO6.45 32 the pres- 
ence of an underlying quantum critical point (QCP) has been conjectured in the doping range 
6 ~ & - 10%. 33 ' 34 This conjecture could be tested in a rather direct way using Raman scattering 
in the superconducting state, see Figs. [7J and [9j The measurement of the Nernst coefficient by 
Daou et al. 36 determined the doping dependence of the JFSD instability in the region of 11 - 18 
% doping and suggested that the pseudogap temperature T* corresponds to the onset of the JFSD 
instability. Raman scattering can directly measure JFSD fluctuations and instabilities generated 
by them, see Figs. [6] and [81 and thus prove the consistency of transport and light scattering data. 
Moreover, the resistivity measurement by Daou et al. 61 suggested that T* goes down to zero in 
the overdoped region, implying the presence of a QCP associated with the JFSD instability inside 
the superconducting state. Theory 3 - 3 -* 3 ^- and transport measurements however, conjecture quite 
different values for the position of the QCP as a function of doping which also could be clarified 
by Raman scattering in the superconducting state, see Figs. [7J and [9j The neutron scattering exper- 
iments for YBC06,45, 32 YBC06.6, 29 ' 30 and YBCO^.st^- suggested a delicate interplay between the 
tendency towards a JFSD and the singlet pairing formation in agreement with theory. 31 ' 35 Raman 
scattering around T c or the pseudogap temperature T* can directly reveal how the JFSD competes 
with the singlet pairing at different doping levels. Available Raman scattering data 62 for YBCO 
with 10 % carrier doping do not suggest the strong enhancement of the low energy spectral weight, 
seen in Fig. [6l but the data were obtained only at a few temperatures. More detailed experimental 
studies including doping dependence are worth performing. 

The La-based cuprate superconductors were extensively discussed in terms of the charge-stripe 
order. 63 However, the scenario based on the JFSD instability was also proposed . 1164 ' 65 Although the 

interpreted the Bi g Raman scattering spectra for La-based cuprates with 



authors of Refs. |62j and 
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10 % in terms of charge stripes, their data exhibit a spectrum very similar to Fig . [(J a), 



62 



and 



indicating 



66 are worth 



direct evidence of the development of JFSD correlations. The data in Refs. 
reconsidering. 

The JFSD is a generic feature in correlated electron systems and occurs in the t-J-— and 
Hubbar d 4 ' 5 ' 67 ' 68 models, in systems where electrons interact via a central forc o 69 ' 70 , and quite gen- 
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erally in Fermi liquids with a van Hove saddle point. Therefore the JFSD instability can be 
expected to occur in a variety of materials. In order to apply the present theory in the normal state 
we had to include self-energy effects. While quantitative features of the Raman spectrum certainly 
depend on details of the self-energy, it is not unreasonable to assume that its qualitative features 
associated with the proximity of the JFSD instability are rather robust. In this sense we hope that 
our results will serve to analyze Raman scattering data for various materials, which possibly lie 
close to a JFSD instability. In particular, compelling but indirect evidence for a JFSD instability 
has accumulated in Sr327 both experimentally^— and theoretically.— - — It would be desirable to 
perform also Raman scattering measurements in this system to confirm the JFSD instability in a 
more direct and decisive way. 

In summary, we have studied Raman scattering from electrons and phonons in the normal 
and superconducting state near a <iFSD instability. In the normal state the inclusion of the elec- 
tronic self-energy is vital for which we have used experimental input from ARPES data in high-T,. 
cuprates. 39,40 Approaching the JFSD instability from the normal state a central peak emerges both 
in electronic scattering and in the spectral function of a phonon with B\ g symmetry. Approach- 
ing the JFSD instability in the superconducting state by decreasing doping concentrations a sharp 
soft mode appears in electronic Raman scattering. This soft mode also appears in the spectral 
function of the phonon as an additional low-energy feature whereas the usual phonon peak is 
nearly unaffected by the proximity of the instability. Our study was motivated by recent trans- 
port measurements 36 which suggest a JFSD instability in a wide doping region in YBCO. Since 
Raman scattering measures directly the correlation function of order parameter fluctuations asso- 
ciated with the JFSD instability such measurements, together with our theoretical curves, could 
confirm in a rather direct way the JFSD in real systems. 
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